Finite size scaling of the correlation length above the upper critical dimension 
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We show numerically that correlation length at the critical point in the five-dimensional Ising 
model varies with system size L as L^''*, rather than proportional to L as in standard finite size 
scaling (FSS) theory. Our results confirm a hypothesis that FSS expressions in dimension d greater 
than the upper critical dimension of 4 should have L replaced by L''^* for cubic samples with periodic 
boundary conditions. We also investigate numerically the logarithmic corrections to FSS in d = 4. 
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I. INTRODUCTION 

Finite size scaling^-^^ (FSS) has been extremely use- 
ful in extrapolating numerical results on finite systems 
in the vicinity of a critical point to the thermodynamic 
limit, in order to get information on critical singularities. 
The basic hypothesis of FSS is that the linear size of the 
system L enters in the ratio i/^oo where S,oo is the cor- 
relation length of the infinite system (which we will call 
the "bulk" correlation length for convenience) and which 
diverges as the critical temperature is approached like 



where 



t = 



Cot' 



T-Tr. 
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measures the deviation from criticality. Here cq is a non 
universal "metric factor''^ and we use the symbol ~ to 
signify "asymptotically equal to" . Hence, if a quantity 
X diverges in the bulk like the FSS form for the 

behavior of X is 



Xo VCoo 



(3) 



where X^ and ci are non- universal scale factors, and ± 
refers to t ^ 0. The scaling fmictions and X are uni- 
versal . In the last expression in Eq. ||2J) we have taken 
the argument of the function in the first expression 
to the power l/v, in order that temperature appears lin- 
early. This has the advantage that a single smooth func- 
tion X, applies both above and below Tc, whereas two 
functions P^ are needed in the first expression in Eq. 

It is often convenient to consider dimensionless quan- 
tities, because these are expected to have Ux = Q- Two 
commonly studied examples are (i) the "Binder ratio"—. 
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where m is the order parameter, and (ii) the ratio of the 
correlation length of the finite system S.L to the system 
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The definition of is not unique (though any reasonable 
definition will give the same scaling form). We shall give 
one definition, which is often used in numerical work, in 
the next section. Again, the scaling functions, W^,g, 
and ^, are universal. 

Note from Eqs. I^J and ||SJ) that, for dimensionless 
quantities like g and S,l/L, data for different sizes in- 
tersect at the critical temperature. Hence dimensionless 
quantities are very convenient because they locate the 
critical temperature in a simple way, from the crossing 
point, without needing to know the values of other quan- 
tities such as exponents. Furthermore, since the scaling 
functions 'g{x) and S^{x) are universal the values of g and 
^l/L at the crossing point (i.e. at Tc) are also universal. 

Finite size scaling, as represented here by Eqs. 
is expected to be valid in the limit L — > oo,t — > 0, with 
L^/'^t arbitrary. Originally proposed on phenomenologi- 
cal grounds, a justification for FSS was later provided by 
Brezini using renormalization group (RG) arguments, at 
least for the case of systems without disorder (which is 
the only case we discuss here). However, Brczin^ also 
noted that FSS breaks down at the "upper critical di- 
mension" du = 4. For d > 4 the critical exponents are 
given by mean field theory, e.g. u = 1/2, and the cor- 
responding field theory is a free theory (i.e. the fluctua- 
tions are Gaussian) since the effective coupling constant 
vanishes at long length scales. This coupling constant is 
irrelevant in the RG sense, but singularities occur when 
it tends to zero, and so it cannot simply be set to its 
"fixed point" value of zero. It is the singularities which 
come from this dangerous irrelevant variable that lead to 
a breakdown of the FSS expressions in Eqs. for 
d > 4. 

Nonetheless, since the bulk behavior for o? > 4 is triv- 
ial, one might imagine that, in this limit, the size depen- 
dence can be also be expressed in fairly simple way and 
this turns out to be the case. As seems to be at least 
implicit in much of the earlier workSiSiSiifliii we make the 
hypothesis that, for cubic samples with periodic bound- 
ary conditions: 
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for d > 4, FSS formulae can still be applied 
but with the system size L replaced by a 
larger length^^ 

where Ai is non-universal. 



(6) 



We shall see that physically i is the correlation length at 
the critical point. With this replacement, Eqs. 
become (remember d > A here) 



(7) 



W^M^) ^~g(c,L'^/H), (8) 



i (c^L^'H) (9) 



with C2 non- universal, where we have noted that v = 1/2 
for d > 4. As before, the scaling functions are universal, 
so the value of g at the crossing point at Tc is universal. 
Furthermore, this universal value has been calculated2i^. 
We see that, at criticality, is of order L'^^'^ which is 
much greater than L for large sizes, a result which, at 
first, seems surprising. The value of S,l/L'^^'^ at critical- 
ity, however, is non universal because of the factor of Ai 
in Eq. This factor occurs because I has dimensions 
of length, and so, for Eq. © to be dimensionally correct, 
Ai must be proportional to a^~'^^^, where a is a micro- 
scopic length scale, e.g. the lattice spacing. Quantities 
involving microscopic length scales are not universal and 
so Ai is not universal. 

There has been extensive discussior&ii*i2ii^iiSii& as to 

whether Eq. © applies to the five-dimensional Ising 
model in the limits L ^ oo, i — > 0. Apparently it 
does^^^, though there appear to be several corrections 
to FSS which conspire to give a "crossing" for small sizes 
at a value of g which differs from the calculated universal^ 
value. 

As noted above, a surprising feature of Eq. lO is that 
the correlation length of the finite system at the critical 
point is greater than the system size. To our knowledge 
there does not appear to have been any direct verifica- 
tion of this prediction for c? > 4 by numerical simula- 
tions. In this paper, we confirm the prediction in Eq. 
by Monte Carlo simulations on the fivc-dimensional Ising 
model. We also carry out similar simulations for the four- 
dimensional Ising model, for which logarithmic correc- 
tions to standard FSS are expected^. 

In Sec. ^ describe the model and some aspects of the 
simulations. The results in five dimensions are presented 
in Sec. mil and the results in four dimension are presented 
in Sec. IIVI We summarize our results in Sec. 
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FIG. 1; Data for S,l/L in d = 5. Clearly the data do not 
intersect at a common point, as would be expected if the 
conventional FSS expression, Eq. applied. 

II. THE MODEL 



The Hamiltonian is given by 

n 



J ^ SiSj , 



(10) 



where the Ising spins take values Si — ±1, and the sites 
i are on a hypercubic lattice in d dimensions of size N — 
L'^. We take d = 4 and 5, and apply periodic boundary 
conditions. The sum is over nearest neighbor pairs of 
sites, and from now on we set J — 1. 



The magnetization per spin is given by 



1 ^ 



(11) 



and the Binder ratio is then given in terms of moments of 
m by Eq. Q . The correlation length of the finite system, 
is given by the following finite difference expression^ 



1 



/ C(kn.in) ^ 



2sin(fc,„in/2) V C(0) 



(12) 



where 
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FIG. 2: Data for ^l/L^^^ in d = 5. Clearly the data intersect 
close to a common point, as expected for the modified FSS 
expression in Eq. Q. The vertical line is at T = 8.7785 which 
is our best estimate for Tc. 



is the Fourier transform of the spin-spin correlation func- 
tion, and kniin — (27r/L)(l, 0, 0) is the smallest non-zero 
wave vector on the lattice. Above and for L oo, 
Eq. H12|l gives the usual second moment definition of the 
correlation length. 

We perform Monte Carlo simulations using the Wolff''^ 
cluster algorithm to reduce the effects of critical slowing 
down. 



III. RESULTS IN 5 DIMENSIONS 
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FIG. 3: A scaling plot for the data in Fig. |21 according to the 
second expression in Eq. Q with Tc = 8.7785. 



presumably because of corrections to FSS, the data scales 
well with Tc = 8.7785. By considering different choices 
for Tc we estimate that Tc = 8.7785(5), consistent with 
the more accurate result 8.77844(2) in Ref. p^ . 



For completeness we also show results for the Binder 
ratio in Fig. 2| As found in other worli^4^4'^4^4^, the 
data for small sizes intersect at a value of g larger than 
the predicted^'^ universal value of 0.4058 • • •. The data 
for larger sizes have intersections at somewhat smaller 
values and presumablj*i& would reach the universal value 
for L oo. 



Data for /L is shown in Fig. ^ for sizes 4 < L < 
16. According to standard FSS, Eq. lO, the data would 
intersect at a common point which is clearly not the case. 
However, according to the modified FSS expression in 
Eq. (PJ, it is data for S,l/L^^'^ which should intersect at 
a common point, and Fig.|21shows that this works pretty 
well. Fig. 121 therefore provides convincing evidence that 
the correlation length at the critical point varies as L^/** 
in 5 dimensions, rather than being proportional to L as 
would be expected in standard FSS. 

A scaling plot of the data in Fig. El according to the 
second expression in Eq. ^ is shown in Fig. |3| Note 
that in addition to the vertical axis being scaled by L^/^, 
rather than L as in standard FSS, T — Tc is scaled by 
i^/^ rather than L^/" (= L^) as in standard FSS. Apart 
from L = 4, for which the data is consistently too low 



IV. RESULTS IN 4 DIMENSIONS 



In four dimensions, Brezini argued that oc 
L{\ogLy/'^ at criticality, and so we expect that FSS ex- 
pressions should be modified by the replacement 

£ = A2L{lnLy/'^ {d ^ i) . (14) 

In Fig. El we show a plot for £,l/L (i.e. without the loga- 
rithmic factor). Clearly the data does not show a com- 
mon intersection. However, including the logarithmic 
factor, the plot in Fig. shows a good intersection with 
only small corrections to FSS. The factor InL can be re- 
placed by ln(L/Lo) where Lq is a microscopic scale, and 
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FIG. 4: Data for the Binder ratio in d = 5. The vertical 
dashed hne corresponds to T = 8.7785 which is our best es- 
timate of Tc from the correlation length data, see Fig. I^. 
The horizontal dashed line corresponds to g = 0.4058 • • •, the 
predicted^*^ universal value. 
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FIG. 5: Data for ^l/L in d = 4. According to conventional 
FSS, Eq. ||KJ, the data should have a common intersection. 
This is clearly not the case. 



FIG. 6: Data for ^L/(i In L)^/* in d = 4. The data intersect 
at close to a common point. 

with an appropriate choice of Lo we get sharper intersec- 
tions. However, ln(L/Lo) = (lnL)(l + InLo/lnij) and 
so including Lq corresponds to an additive correction to 
FSS (which vanishes only logarithmically). It is difficult 
to separate this from other corrections to FSS, and so we 
don't feel we can give a reliable estimate for Lq- 



V. CONCLUSIONS 

We have demonstrated that the FSS behavior of the 
correlation length (for a cubic sample with periodic 
boundary conditions) in five dimensions follows Eq. 
which is the expected modification of FSS for the case 
d > A. This provides confirmation that the standard 
FSS expressions, e.g. Eqs. ©-(O, can be simply mod- 
ified above c? = 4 by the replacement^^ L — > £ oc L''/*, 
which gives Eqs. Q-®. This had been verified before 
for the Binder ratio, but not, to our knowledge, for the 
correlation length. It is interesting that the correlation 
length at the critical point is of order i and hence much 
bigger than the system size L. This is possible because 
the long wavelength fluctuations arc non-interacting near 
criticality for d > 4. We also demonstrated the expected 
logarithmic modification to FSS of the correlation length 
for d precisely equal to 4. 

It is also interesting to ask what are the corresponding 
results with d > 4 for other geometries and boundary 
conditions. For the "strip" geometry, where the sample 
is infinite in one direction and of size L in the others, 
Brezin^ showed that the correlation length at the critical 
point varies as L^'*"^)/'^ (which is reasonable since FSS is 



5 



done only with respect to the d—1 finite dimensions). It 
is then natural to expect that FSS will then work with L 
replaced throughout by L(''~i)/3_ 

For free boundary conditions, it seems obvious that 
even for d > 4 the behavior of the system will be af- 
fected when becomes of order L, rather than only 
change when becomes of order the much larger length 
£. Hence we expect that the standard FSS expressions, 
Eqs. ©-HI would apply with v = 1/2. The ratio ^oc/£ 
may also enter but, since £ ^ L for large L, such terms 
would presumably be corrections to the scaling terms 
which involve ^oo/L. Since FSS for models with free 



boundary conditions in d > 4 is poorly understood, it 
would be interesting to investigate such models in some 
detail. 
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